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Abstract 

We present a robust method for quantum process tomography, which yields a set of Lindblad 
operators that optimally fit the measured density operators at a sequence of time points. The 
benefits of this method are illustrated using a set of liquid-state Nuclear Magnètic Resonance 
(NMR) measurements on a molecule containing two coupled hydrogen nuclei, which are sufficient 
to fully determine its relaxation superoperator. It was found that the complete positivity constraint, 
which is necessary for the existence of the Lindblad operators, was also essential for obtaining a 
robust fit to the measurements. The general approach taken here promises to be broadly useful 
in studying dissipative quantum processes in many of the diverse experimental systems currently 
being developed for quantum information processing purposes. 



PACS numbers: 03.65.Wj, 03.65.Yz, 07.05.Kf, 33.25. +k 
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INTRODUCTION 



An important task in designing and building devices capable of quantum information 
processing (QIP) is to determine the superoperators that describe the evolution of their 
component subsystems from experimental measurements. This task is commonly known in 
QIP as Quantum Process Tomography (QPT) []]]. The superoperators obtained from QPT 
allow one to identify the dominant sources of decoherence and focus development efforts 
on eliminating them, while precise knowledge of relevant parameters can be used to design 
quantum error correcting codes the and/or decoherence- free subsystems that circumvent 
their effects [[I], H . Methods have previously been described by which the "superpropagator" 
V of a quantum process can be determined || ||. Assuming that the process' statistics are 
stationary and Markovian |5|, |6| , a more complete description may be obtained by determin- 
ing the corresponding "supergenerator" , that is, a time-independent superoperator Q from 
which the superpropagator at any time t is obtained by solving the differential equation 
V(t) = —ÇV(t). The formal solution to this equation is exp(— Qt), where "exp" is the 
operator exponential. 

The purpose of this paper is to describe a data fitting procedure by which estimates of a 
supergenerator can be obtained. This problem is nontrivial because, as in many other data 
fitting problems, the estimates obtained from straightforward (e.g. least-squares) fits turn 
out to be very sensitive to small, and even random, errors in the measured data. In some 
cases, this may result in a supergenerator that is obviously physically impossible; in others, 
it may simply result in large errors in the generator despite it yielding a reasonably good 
fit to the data. Parameter estimation problems with this property are commonly known 
as ill-conditioned 0, |§. The main result of this paper is that, although the problem of 
estimating a supergenerator from measurements of the superpropagators at various times is 
ill-conditioned, this ill-conditioning can be greatly alleviated by incorporating prior knowl- 
edge of the solution into the fitting procedure as a constraint. The prior knowledge that we 
use here is a very general property of open quantum system dynamics known as complete 
positivity. 

Roughly speaking, complete positivity means that if V is a superoperator that maps a 
density operator of a system to another density operator, then any extension of the form 
X (gi V also returns a positive operator, where X denotes the identity map on the extension 
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of the domain. The most general form of a completely-positive Markovian master equation 
for the density operator p of a quantum system is known as the Lindblad equation || ||. 
This may be written as 

M 

pit) = - t [H, p(t)} + \ {[L m p(t), Li] + [L m , p(t)Ll\) , (1) 

m=l 

where h — 1, t is time, H is the system's Hamiltonian, the L m are known as Lindblad 
operators, and the Ü m denote their adjoints. It is easily seen that the Lindblad equation 
necessarily preserves the trace tr(p) = 1 of the density matrix, meaning tr(p) = 0, and a 
little harder to show that it also preserves the positive-semidefinite character of the density 
operator. Proofs that it has the yet-stronger property of complete positivity may be found 
in Refs. H |9], [TÜ], [TT}| . The QPT method we describe in this paper relies upon the Lindblad 
characterization of complete positivity. 

The paper is organized as follows. In the first part of the paper we present a computa- 
tional procedure which fits a completely positive supergenerator to a sequence of estimates of 
the superpropagators of a quantum process at múltiple time points. This procedure initially 
extracts an estimate of the decoherent part of the supergenerator, without the Hamilto- 
nian commutation superoperator (which is assumed to be available from independent prior 
measurements). It then refines this initial estimate via a nonlinear least-squares fit to the 
superpropagators, in which complete positivity is enforced by adding a suitable penalty func- 
tion to the sum of squares minimized. Finally, any residual non-completely-positive part 
of the supergenerator is "filtered" out by a matrix projection technique based on principle 
component analysis flU , |P2 . 



In the second part of the paper, the procedure is validated by using it to determine the 
natural spin-relaxation superoperator of a molecule containing two coupled spin 1/2 nuclei in 
the liquid state from a temporal sequence of density operators. These in turn were derived by 
state tomography, meaning a set of NMR measurements sufficient to determine the density 
matrix. In the process we confirm the ill-conditioned nature of the problem, and that the 
complete positivity constraint is needed to obtain a robust estimate of the supergenerator. 
The final results are used to compute the corresponding Lindblad operators, but these were 
difficult to intèrpret. Hence the Hadamard representation of T2-relaxation dynamics 0] was 
used to derive a new set of Lindblad operators which are easier to intèrpret and explain most 
of the relaxation dynamics. The information these operators convey agree with theoretical 



expectations as well as with some additional independent measurements, in support of the 
accuracy of the results obtained. 



COMPUTATIONAL PROCEDURE 

In this paper we are concerned with an iV-dimensional open quantum system (N < oo), 
the dynamics of which are described by a Markovian master equation of the form |5|, |6|, [Ï4j : 

% = -1-[Hl P) ~ ^(P-Peq) 

dí (2) 

^ d ^ A n-i <r> 

—rr = -ïnpA~i<-PA 
dt 

In this equation, h = 1, p — p(t) is the system's density operator, p eq this operator at 
thermal equilibrium, p& = p — p eq , H is the system's internal Hamiltonian, TC the corre- 
sponding commutation superoperator, and TZ is the so-called relaxation superoperator. The 
equivalence of the first and second lines follows from the fact that p eq is time-independent 
and proportional to the Boltzman operator exp(—H/kQT), so that it commutes with H. 
By choosing a basis for the "Liouville space" of density operators, the equation may be 



written in matrix form as |TT], [14], [15 

^ = -(*2Í + £)Ua> = -QUa) (3) 

where the underlines denote the corresponding matrices and \_p) is the iV 2 -dimensional col- 
umn vector obtained by stacking the columns of the density matrix _p on top of each other 
in left-to-right order||TT|. A numerical solution to this equation at any point t in time may 
be obtained by applying the propagator Vit) to the initial condition \j)/^{0)), where the 
propagator is obtained by computing the matrix exponential V_(t) = exp ( — Ç_t) [ 15] , [Tü . 
Note that 2Í and ]Z do not commute in general, and that the sum Ç_ = iTí + TZ will not 
usually be a normal matrix (one which commutes with its adjoint). This in turn greatly 



reduces the efhciency and stability with which its matrix exponential can be computed (T^ 



(although this was not an issue in the small problems dealt with here), and we expect it to 
also significantly complicate the inverse problem. 

In this section we describe an algorithm for solving this inverse problem, that is to 
determine the relaxation superoperator TZ from an estimate of the Hamiltonian H_ together 



with estimates of the propagator V m = V_(t m ) at one or more time points ti < t 2 < ■ ■ ■ < ím- 
This problem, like many other inverse problems, turns out to be ill-conditioned, meaning that 
small experimental errors in the estimates of the H_ and V m may be amplified to surprisingly 
large, and generally nonphysical, errors in the resulting superoperator TZ R]. For example, 
if one tries to estimate TZ in the obvious way as 

U « (-iU-\o K {V l ))/t l , (4) 

one will generally obtain nonsense even if Jt and V\ are known to machine precision, be- 
cause of the well-known ambiguity of the matrix logarithm with respect to the addition of 
independent múltiples of 2 1 tc onto its eigenvalues. Using the principal branch of the loga- 
rithms will only work if 7í is small compared to TZ, and the only reasonably reliable means 
of resolving the ambiguities is to utilize additional data and/or prior knowledge of the so- 
lution. Even then, a combinatorial search for the right múltiples of 2 1 tt may be infeasibly 
time-consuming. 

An alternative to the logarithm which utilizes data at múltiple time points and is capable 
of resolving the ambiguities even when Ti is much larger than TZ is to estimate the derivative 
at t — of 

e ü ^ 2 V(t)e ü ^ 2 = e"*- + 0(t 2 ) . (5) 

This derivative is obtained by Richardson extrapolation using central differencing about 
t = over a sequence of time points such that t m = 2 m ~ 1 ti (m = 1, . . . , M), according to 
the well-known procedure ||18|| : 



for m from 1 to M do 

T) — om-2 ( p it m K]2 -p p it m ni2 _ p -it m H/2 <r>-l „-it m H/2\ . 

±±l,l+M-m ■— * [L·. L·mL· L·. L·m L·. ) ) 

for £ from 1 to m — 1 do 

Ql+l,l+M-m := ^-(.,1+M-m + {Qe,l+M-m ~ ^-i,M-m) I ~ 1) i 

end do 
end do 

The inverse = V_{—t m ) is assured of existing unless long times are used or the errors 
in the data are large. The method produces an estimate of the derivative at t — that 
is accurate up to 0(tl M ), and which may be increased by computing the exponential from 
the highest-order estimate at further interval halvings. The method performs well when the 



relative errors in the Hamiltonian 5 <C I/^AvÍm), where Au is the range of freqüències 
present in the Hamiltonian, but it tends to emphasize the errors in V l rather than averaging 
over the errors at all the time points. Hence we do not recommend that it be used alone, 
but rather as a means of obtaining a good starting point for a nonlinear fit to the data, as 
will now be described. 

This nonlinear fit involves minimization of the sum of squares 

A II 2 

x 2 {mu,v u ...,v M ) = J2 \\^(-^2í+mtm) -£ m f (6) 



m=l 



with respect to 71, where || ■ ||p denotes the squared Frobenius norm (sum of squares of 
the entries of its matrix argument). Previous results with similar minimization problems 
indicate that \ 2 will have many local minima []TÜ| , making a good starting point absolutely 
necessary (even disregarding the 2zn ambiguities discussed above). The derivatives of this 
function may be obtained via the techniques described in Ref. but the improvements 
in efficiency to be obtained by their use are likely to be of limited value in practice given 
that all the resources needed, both experimental and computational, grow rapidly with TV 
(which itself grows exponentially with the number of qubits used in quantum information 
processing problems). In addition, the quality of the results matters a great deal more than 
the speed with which they are obtained, and the quality will not generally depend greatly 
on the accuracy with which the minimum is located. 

For these reasons, we have used the Nelder-Mead simplex algorithm ||20|| , as implemented 
in the MATLAB™ program, for the small (two qubit) experimental test problem described 
in the following section. This has the further advantage of being able to avoid local minima 
better than most gradient-based optimization algorithms. Preliminary numerical studies, 
however, exhibited the anticipated ill-conditioning with respect to small perturbations in 
the data, even when TZ was constrained to be symmetric (implying a unital system which 
satisfies detailed balance) and positive semidefinite (as required for the existence of a finite 
fixed point). Therefore it is necessary to incorporate additional prior information regarding 
TZ into the minimization. The information that we have found to be effective is a property 
of TZ known as complete positivity [|5|, |6|, |ÏT . 



An intrinsic definition which does not involve an environment was first given by Kraus 



21 , and states that a superpropagator V is completely positive if and only if it can be 
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written clS cL "Kraus operator sum" , namely 

N 2 

V(p) = Y, K ïP K \, (7) 

e=i 

where p = p^ and K\, . . . , K^i all act on the system alone. Another intrinsic definition 
subsequently given by Choi states that a superpropagator is completely positive if and 
only if, relative to any basis of the system's Hilbert space, the so-called Choi matrix is 
positive semidefinite [nj|, namely 

N-l N-l 
i,j=0 i,j=0 

This equation uses the notation of quantum information processing, in which the i-th ele- 
ment ary unit vector is denoted by \i) (0 < i < N), V( \ i)(i \ ) is the N x N matrix of the 
operator obtained by applying the superpropagator V to the projection operator given by 
\i){j\ versus our choice of basis, and V_ is the N 2 x iV 2 matrix for V versus the Liouville 
space basis \i)\j) (as for Jt etc. above). It can further be shown that the eigenvectors kg 
of the Choi matrix are related to the (matrices of an equivalent set of) Kraus operators by 
\K e ) = ^/k~£ kg, where Kg > are the corresponding eigenvalues and the "ket" \Kg) indicates 
the column vector obtained by stacking the columns of Kg on top of one another JTTJ . 

These results can be used not only to compute a Kraus operator sum from any completely 
positive superpropagator given as a "supermatrix" acting on the iV 2 -dimensional Liouville 
space, but also to "filter" such a supermatrix so as to obtain the supermatrix of the com- 
pletely positive superpropagator nearest to it, in the sense of minimizing the Frobenius norm 



of their difference This is done simply by setting any negative eigenvalues of the Choi 



matrix to zero, rebuilding it from the remaining eigenvalues and vectors, and converting this 
reconstructed Choi matrix back to the corresponding supermatrix. Although this generally 
has a beneficial effect upon the least-squares fits versus x 2 ( as defined above), it is still 
entirely possible that the sequence of filtered propagators V' m will not correspond to a com- 
pletely positive Markovian process, so that no time-independent relaxation superoperator 
1Z can fit it precisely. This, together with the ill-conditioned nature of the problem, implies 
one will still not usually obtain satisfactory results even after filtering. For this reason we 
shall now describe how the above characterizations of completely positive superpropagators 
can be extended to supergenerators. 



As indicated in the Introduction, completely positive Markovian processes, or quantum 
dynamical semigroups as they are also known, may be characterized as those with a generator 
Q that can be written in Lindblad form [EL |6], |9], |Ï0|, |ÏTJ. On expanding the commutators in 
Eq ([!]), this becomes 

N 2 

-Q{P) = ~t[H, p] + \ ( 2L mP L l ~ L L· L mP ~ P L L· L m) ■ (9) 

m=l 

The operators L m are usually called Lindblads. It may be seen that the Choi matrix C 
associated with —TZ is never positive semidefinite, because (L\Q\L) — tij—TZ) < 0. Never- 
theless, it can be shown that any trace-preserving TZ (meaning tr(7£(£)) = 0) is completely 
positive if and only if a certain projection of Ç is positive semidefinite, namely £_€_£_ where 



£_ = L ® L ~ \L) (L\ [ 11 1 . In this case an equivalent system of orthogonal Lindblads is 
determined by \L m ) = y/X^í m , where A m > are the eigenvalues and í m the eigenvec- 
tors of S_ÇS_. In the event that £_€_£_ has negative eigenvalues we can simply set them to 
zero to obtain a similar but completely positive supergenerator, much as we did with the 
Kraus operators. Most importantly, however, this characterization of completely positive 
supergenerators gives us a means of enforcing complete positivity during nonlinear fits to a 
sequence of propagators at múltiple time points. 

The following section describes our experience with applying this approach to a sequence 
of propagators obtained from liquid-state NMR data. The complete positivity of the re- 
laxation superoperator was maintained by adding a simple penalty function onto the sum 
of squares that was minimized by the simplex algorithm, as described above. This penalty 
function consisted of the sum of the squares of the negative eigenvalues of the corresponding 
projected Choi matrix. While more rigorous and efficient methods of enforcing the projected 
Choi matrix to be positive semidefinite are certainly possible, this strategy was suflicient to 
demonstrate the main result of this paper, which is that the complete positivity constraint 
greatly alleviates the ill-conditioned nature of such fits. 



EXPERIMENTAL VÀLID ATION 

The experiments were carried out on a two-spin | system consisting of the hydrogen 
àtoms in 2,3-dibromothiophene (see Fig. 1) at 300K dissolved in deuterated acetone, using 
a Bruker Avance 300 MHz spectrometer. The internal Hamiltonian of this system in a frame 
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FIG. 1: Molecule of 2,3-dibromothiophene with the two protons labeled 1 and 2. The chemical 
bonds among the àtoms are indicated by double parallel lines, and a transparent "dot-surface" 
used to indicate their van der Waals radii. 

rotating at the frequency of the second spin is 

H = H mt = v{v x o\ + f a 1 ■ a 2 ) (10) 

where V\ = 161.63Hz is the chemical shift of the first spin, J = 5.77Hz is the coupling 
between the spins, and cr = [<j x , a y , a z ] are Pauli spin operators. 

The "quantum operation" we characterized was just free-evolution of the system un- 
der its infernal Hamiltonian, together with decoherence and relaxation back towards the 
equilibrium state p eq ~ o\ + a 2 . In liquid-state NMR on small orgànic molècules like dibro- 
mothiophene, this process is mediated primarily by fluctuating dipolar interactions between 
the two protons as well as with spins neighboring molècules, and since the correlation time 
for small molècules in room temperature liquids is on the order of picoseconds, the Marko- 
vian approximation is certainly vàlid JÏ4], We add that our sample was not degassed so 
that the presence of dissolved paramagnetic O2 shortened the T\ and T 2 relaxation times. 

The experiment consisted of preparing a complete set of orthogonal input states (that 
is, density matrices), letting each evolve freely for a given time T, and then determining 



the full output states via quantum state tomography |25], ^6|. Since only "single quantum" 
coherences can be directly observed in NMR [pü| , ^M , this involves repeating the experiment 
several times followed by a different readout pulse sequence each time, until all the entries of 
the density matrix have been mapped into observable ones. The experiments were carried out 
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index operator order 
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4 
5 
6 
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4x4 identity 
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a x aç - a y a x 
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10 
11 
12 
13 
11 
15 
16 



1 2 

a y a z 



1 2 

a z a y 
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(J V (T. 
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TABLE I: Table of operators (versus Cartesian basis) in the transition basis used for the density 
and superoperator matrices, the corresponding matrix indices and their coherence orders (see text). 

at four exponentially-spaced times T, as required by the Richardson extrapolation procedure 
described above, specifically T = 0.4, 0.8, 1.6 and 3.2s. 

To describe the density and superoperator matrices, the so-called "transition basis" was 
used ||14|| . This Liouville space basis is intermediate between the Cartesian basis and the 
Zeeman (or polarization and shift operator JL4J]) basis, in that the basis elements are all 
Hermitian like those of the Cartesian basis, but like the Zeeman basis they have a well- 
defined coherence order, or difference in total angular momentum along the applied magnètic 
field Bq between the two Zeeman states connected by the transition. These basis states are 
listed in TABLE I versus the Cartesian basis. 

This basis was chosen because the relaxation superoperator 1Z is expected to have the 



"Redfield kite" structure in this basis jTJJ]. This block diagonal structure arises because 
the difference in frequency between transitions of different coherence orders, given that 
the Zeeman interaction dominates all others, is large enough to average out these other 
interactions including those responsible for decoherence and relaxation, in effect decoupling 
the blocks from one another so that no cross relaxation occurs between them (see Fig. 2). 
This so-called "secular approximation" considerably reduces the number of parameters in 
the superoperator from 256 = ((2 2 ) 2 ) 2 to 81 = 3 2 + 2 2 + 8 2 + 2 2 (since neither the identity 
nor the other diagonal (cr z ) basis elements are expected to cross relax with any non diagonal 
elements). 
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Population 
biock (4x 4) 

Zero quantum 

block (2x 2) 



Singla quartum 
blook (8x8) 



Double quantum 

block (2x 2) 

FIG. 2: Redfield kite struture of the relaxation superoperator expressed in the transition basis. 
The non-shaded area corresponds to Blocks of different coherence order are decoupled from each 
other. 

An additional reduction may be obtained by assuming detailed balance: the microscopic 
reversibility of all cross relaxation processes. The relaxation superoperator reconstructed 
from the experimental data was bordered with an initial row and column of zeros to force 
7Z(I) = 0, because the totally random density matrix J/4 cannot be observed by NMR 
spectroscopy and because of the trace-preserving process assumption that we made. This 
may be done providing TZ operates on pa = P — p eq , and together with detailed balance it 
implies that the supermatrix TZ will be symmetric, reducing the number of parameters to 
be estimated to only 48 = 6 + 3 + 36 + 3. 

The result of applying the fitting procedure without the complete positivity constraint 
to the initial estimate obtained by Richardson extrapolation is shown in FIG. 0(a). It may 
be seen that the rates did not vary in a systematic fashion with the coherence order, and 
that large cross-relaxation rates were found, which is not consistent with the physics of 
spin relaxation in liquid-state NMR spectroscopy. In addition, this relaxation superoperator 
implies that spin 1 has a T 2 ~ 2.3s while spin 2 has a T 2 ~ 4.6s, in disagreement with the 
independent measurements of T 2 given below. 

The fit after adding the complete positivity constraint is shown in FIG. |3|(b), again start- 
ing from the results of Richardson extrapolation. It may now be seen that the results do 
vary systematically with coherence number, and that the resulting relaxation superoperator 
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FIG. 3: Three different estimates of the relaxation superoperator of 2, 3-dibromothiophene in the 
transition basis, indexed as indicated in Table I. (a) Relaxation superoperator obtained from a 
least-squares fit, without the complete positivity constraint, of the exponential exx> (—z(TÍ+TZ)t m ) to 
the propagators V m at the corresponding times (íi = 0.4, Í2 = 0.8, Í3 = 1.6, Í4 = 3.2s) with respect 
to the unknown relaxation superoperator ]Z, starting from the results of Richardson extrapolation 
(see text), (b) The relaxation superoperator obtained from a fit to the same data and with the 
same starting value of ]Z, but with the complete positivity constraint included in the fit. (c) The 
relaxation superoperator obtained by assuming that 7i and 1Z commute, and using the average of 
the estimates obtained by taking the logarithms of the absolute vàlues of the eigenvalues of the 
propagators over all four time points as the final estimate (see text). 

is very nearly diagonal. To obtain further evidence for the validity of this superoperator, 
we measured the single spin 7\ (longitudinal, or a z ) and T 2 (transverse, or o~ x = a y ) re- 



laxation rates, using the weli-established inversion-recovery and CPMG experiments p4 
The results for both spins were T\ = 5.6s and T2 = 2.7s, which agree quite well with the 
vàlues of 5.6s and 2.6s obtained from this relaxation superoperator. Although this is obvi- 
ously a relatively simple relaxation superoperator, it is reasonable to expect that a complete 
positivity constraint will substantially improve the estimates of more complicated super- 
operators containing nonzero cross-relaxation rates that cannot be obtained from Standard 
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experiments. 

Finally, the cross-relaxation rate between the population terms a\ and a^, which is due to 
the well-known nuclear Overhauser effect |27||, is essentially zero in Fig. 3(b). This can occur 
when the overall rotational correlation time of the molecule plus its "solvent-cage" is on the 
order of lns, but was somewhat unexpected given the small size of 2, 3-dibromothiophene. As 
a result, we carried out the selective inversion recovery experiment that consists of inverting 
selectively the longitudinal magnetization of one of the two protons and looking at the 
evolution of the magnetization of the other one while the first relaxes towards thermal 
equilibrium. The change in longitudinal magnetization of the second proton was measured 
to be less than 1% of the unperturbed magnetization essentially revealing no NOE effect 
and providing yet further evidence for the validity of this superoperator. The discrepancy 
between this result and the one given in [^4] is explained by the presence of dissolved 
paramagnetic O2 in our sample so that the Ti relaxation time was shortened in such a way 
that the NOE effect became almost unobservable (our T\ was measured to be 5.6s while the 



T\ reported in [gjj is around 47. 5s). 

Because of the substantial degeneracy of the diagonal elements with the same coher- 
ence order, the superoperator in Fig. 3(b) was also very nearly diagonal in the eigenbasis 
of the Hamiltonian commutation superoperator 7í, so that 7i and 1Z very nearly com- 
mute. This allowed further estimates to be obtained directly from the superpropagators 
P m ~ exp ( — iHtm, ) exp ( —lZt m ) , simply by taking the (real) logarithms of the absolute vàl- 
ues of their eigenvalues and thereby cancelling the phase factors from the Hamiltonian's 
exponential. From FIG. ||](c) we see that the result of averaging these estimates over all four 
evolution times is very similar to the completely positive estimate in FIG. |3](b) (correlation 
coefficient 0.80; squared norms of difference over average 0.90). We note that the estimate in 
FIG. |3](c) did not explicitly assume the Redfield kite structure, thereby providing a further 
consistency check on our results. 



INTERPRETATION VIA LINDBLAD AND HADAMARD OPERATORS 

In this section we present a system of Lindblad operators which act on the density operator 
to give essentially the same derivative as the relaxation superoperator described above (see 
FIG. H). As described in the foregoing " Còmput ational Procedure" section, such a system 
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of Lindblad operators may be obtained by diagonalizing the corresponding projected Choi 
matrix, although it will be seen that a more easily interpreted system was obtained by 
considering the parts of 7Z responsible for T\ and T 2 relaxation separately and by using the 
so-called "Hadamard relaxation matrix" formalism O]. Because these calculations were 
somewhat involved, however, the details of the Hadamard operators calculations are given 
in an appendix. We add that from here on, the relaxation superoperator 1Z will refer to the 
matrix shown in FIG. 0(b). 

These representat ions of relaxation processes are normally applied to the density matrix 
in the Zeeman basis jf£ e (regarded as the computational basis in QIP), which requires 
converting the supergenerator 1Z from the transition to the Zeeman basis. This is easily done 
via a unitary transformation, TV ee = 2 U_]V ra U_ , where U_ \{f ra ) = |£ zee ) (the matrix U_ may 
be derived from TABLE I; the factor of 2 corrects for a change in norm due to the fact that 
the transition basis is Hermitian). Although any relaxation superoperator can be modified 
to act directly on the density operator p rather than its difference with the equilibrium 
density operator pa = P — p eq (vide supra) by taking the right-projection 1Z (X — |p) eq (X|) 
|29| . j3Üf , this makes only a negligible change to 1Z since in liquid-state NMR p eq differs from 
the identity I by < 10 -5 . For that reason, the treatment of T\ relaxation given below was 
considerably simplified by treating it as a unital (identity preserving) process acting on p^. 

As described following Eq. (^), a complete system of Lindblad operators may be obtained 
by diagonalizing the projected Choi matrix 



where it is assumed the eigenvalues have been ordered such that X m > X m +i for m = 
1, . . . , N 2 — 1, and defining the Lindblad matrices such that for all X m > 0: 



This gave rise to a total of 11 Lindblads, the phases of which were chosen so as to make 
them as nearly Hermitian as possible. Once this was done, all 11 were within 2% of being 
Hermitian. 

The relative contributions of these Lindblads to the overall relaxation of the spins can be 
quantified by the squared Frobenius norms ||Iv m || 2 : = A m . This calculation shows that about 
35% of the mean-square noise resided in the first Lindblad, 



£Ç£ = VAV} , 



(11) 




(12) 



0.346 (al + al) + 0.025 o^V 2 , 



(13) 
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which represents strongly correlated dephasing with a T2 for both spins of ~ 4.2s ||13|| , much 
as expected. The next four largest Lindblads together contributed, roughly equally, another 
43% to the total mean square noise, but were considerably more difficult to intèrpret: 



L 2 « -0.013a^ - 0.045aJ - 0.153oJ - 0.061 

+ 0.150 g$c£ - 0.039 ctJctJ + 0.111 ct^ + 0.106a z VJ 

L 3 « + 0.046 - 0.026^ - 0.057 2^ - 0.266 (14) 

L 4 « - 0.024a^ - 0.006aJ - 0.081 - 0.077 

- 0.155 qiq_l - 0.193 aja^ + 0.002 ^ - 0.012 £0% 

L 5 « - 0.017 0^ - 0.060 - 0.071 q_l - 0.090 

+ 0.090 - O.OOIctJ^ - 0.183 0^ - 0.118^ctJ 

It can be shown that L3 contributes about 0.15s _1 to the decay rates of the single-quantum 
coherences (single-spin fiips), bringing down the decay time T 2 w 2.6s and, save for some 
small cross terms among in the single quantum block, rather little else. 

The superoperators corresponding to each of the remaining 9 Lindblads separately all 
contained significant cross terms between the populations and the zero or double quantum 
coherences, in violation of the secular approximation |Ï4]| . Only on summing over all of 
them do these nonphysical cross-terms cancel out, leaving a largely diagonal relaxation 
superoperator behind: the ratio of the mean-square value of the off-diagonal entries of TZ 
to that of the diagonal entries was 1.3% in the transition basis and 3.8% in the Zeeman; 
the latter dropped to 1.8% on excluding the block corresponding to T\ relaxation of the 
populations (vide infra). The nonphysical nature of most of the Lindblads is clearly an 
artifact of the way that our procedure for calculating them forces them to be orthogonal 
and minimal in number. In order to physically intèrpret the dominant relaxation processes, 
we therefore focus our attention first on T\ relaxation among the populations (diagonal 
entries of the density matrix in the Zeeman basis), along with the associated nonadiabatic 
T2 relaxation, and then try to account for the remaining T2 relaxation via simple adiabatic, 
albeit correlated, processes. 

Based on Hadamard operator calculations that can be found in the appendix, and by 
treating the T\ and T 2 relaxation processes separetely, we found that the four Lindblad 
operators which describe the T\ relaxation of the first spin may be replaced by the Hermitian 
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operators: 



= V0.1532 \ al 



V£ = V0.1532|í7j (15) 



= V0.1532±a x V 2 



= VO. 1532 i <x y V 2 

and for the second spin: 

L^ 2 = v / 0.1528i(T 2 , Z,g = V0.1528 i trj, 
L~ = V0.1528 1 (T Z V 2 , L% = v / 0.1528Ía z V y 2 . 



(16) 



In addition, the near-degeneracy of the (1, 4) and (2, 3) rates in the relaxation superoperator 
in the Zeeman basis can be used to combine the associated Lindblad operators into four 
multiple-quantum Ti Lindblad operators based on the average rate: 



r xx 
L· Tl 


= V0.0252 \ o\ol, 


rxy 
L· Ti 


= V0.0252Í 


a x a y 


7-yx 

L· Ti 


= v / 0.0252Ía y V 2 , 


ryy 

L· Ti 


= V0.0252Í 





(17) 



By working through some examples, it may be seen that the sum of the Lindbladian 
superoperators for each of the three sets of four Lindblad operators above also causes all 
the off -diagonal entries of jf£ e to decay with the same rate constant \jT\. This corresponds 
to nonadiabatic T 2 relaxation. Likewise, based on the results above, we substracted the 
nonadiabatic T 2 contribution above to get the adiabatic T 2 contribution and thereby obtained 
the following three Lindblad operators: 



L% = v/09560 ^ K + a 2 ), Lff = V / 0l72Ï ± (a\ - a\) , 



(18) 



L^f = V0Ï93Ï \ o\ o\ . 
These correspond to totally correlated, totally anticorrelated, and pure single- quant um T 2 



relaxation, respectively |ï3fl . As a result, the Hadamard product formalism allowed us to 
obtain a description with a clearer physical interpretation. However, this left a slight dis- 
crepancy between the new results and the original data, because some of the assumptions 
made, such as the near degeneracy of some decay rates mentioned above, were only ap- 
proximations. In order to assess how much of the total experimental superoperator 1Z was 
accounted for by our simplifications, we computed and from the above Lindblad 
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matrices, obtained as described in the Appendix, and computed the relative discrepancy in 
the superoperators, 

||^ zee -DÍM(|^ e ))-Çh oj(l^ e ))ll 2 (ig) 



1 1 1Z zee 1 1 2 

where "zee" indicates that the above quantities were expressed in the Zeeman basis, and 
where Choi is a function which distributes the entries of its 4 x 4 argument over the corre- 
sponding entries of a 16 x 16 matrix and sets all the other entries to zero. This gave a value 
of 6.3%, indicating that the approximations introduced in the Appendix were sufficient to 
account for a large majority of the observed relaxation dynamics. 



CONCLUSION 



In this paper, we have demonstrated a robust procedure by which one can derive a set 
of Lindblad operators which collectively account for a Markovian quantum process, without 
any prior assumptions regarding the nature of the process beyond the physical necessity of 
complete positivity. This procedure should be widely useful in studies of dissipative quantum 
processes. In the appendix, we have further shown how one can use one's physical insight into 
the particular system in question to derive the physical "noise generators" of the system. 
We believe this two-step process is illustrative of how Quantum Process Tomography on 
many physically distinct kinds of quantum processes should be done. 
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APPENDIX 



In this appendix, we shall derive a more compact and direct representat ion of the T 2 relax- 
ation processes operative in 2,3-dibromothiophene via the so-called "Hadamard relaxation 
matrix" JÏ3| (vide infra). 

The populations block of the relaxation superoperator corresponds to indices 1 through 
4 in the transition basis (see Table I.) and the nonzero entries of |X)(Z| i n the Zeeman. The 
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vàlues obtained from the completely positive least-squares fit shown in FIG. ^| are 



iden 



ptra 



0.0000 0.0000 0.0000 0.0000 

0.0000 0.1780 -0.0002 0.0089 

0.0000 -0.0002 0.1784 -0.0097 

0.0000 0.0089 -0.0097 0.3061 



ITT) |U> UT> III) 

0.3301 -0.1435 -0.1617 -0.0249' 
-0.1435 0.3129 -0.0254 -0.1440 
-0.1617 -0.0254 0.3500 -0.1630 
-0.0249 -0.1440 -0.1629 0.3319 . 



pzee 



(20) 



in the transition (left) as well as the Zeeman (right) bases. The matrices R^? and are 
related by the Hadamard transform W_ |T3|, 



W = 



that is 



2WR Z £W 



ptra 

Rt x 



1111 
1 1-1-1 
1-1 1-1 
1-1-1 1 



pzee 



(21) 



(22) 



since W 2 = L In the absence of cross-correlation, symmetry considerations imply that 



R_t^ should be centrosymmetric [p4|, |27|] , and hence we shall use the symmetrized version 
^{BlTi + ^x^x SíTi ílxüx) an d its Hadamard transform in the following, which are 



ptra 



0.0000 -0.0000 -0.0000 0.0000 
0.0000 0.1780 -0.0002 0.0000 
0.0000 -0.0002 0.1784 0.0000 
0.0000 0.0000 0.0000 0.3061 



0.3310 -0.1532 -0.1528 -0.0249 
-0.1532 0.3315 -0.0254 -0.1528 
-0.1528 -0.0254 0.3315 -0.1532 
-0.0249 -0.1528 -0.1532 0.3310 



pzee 
RTx 



(23) 



As noted in the main paper, the T\ of both spins is ~5.6 s, while the NOE rate (connecting 
a\ and a 2 in the transition basis) is negligibly small (—0.0002). 

The entries of —R^ are equal to the diagonal entries of the Choi matrix C of — 7^ zee , and 
the o/f-diagonal entries of R^ e are the only non-negligible entries in their respective rows 
and columns of C_. Therefore, they are eigenvalues of the Choi matrix as well as its projection 
S_ÇS_, and their corresponding eigenvectors are elementary unit vectors \k)\j) relative to the 
Zeeman basis. It follows that the Lindblad operator for the (j, k)-th off-diagonal entry of 
R_t^ may be written as lJ T = (— (j\R^ l^)) 1 ^ 2 1^)01; an d its contribution to ^ e is given 
by 



(j\Rn\k) (\k)(j\2T\j)(k\ - $m\/>? - hr\j)(j\ 



(24) 



The symmetry of R^ implies that the eigenvalues of C corresponding to single spin flips (the 
so-called single-quantum coherences) are four-fold degenerate, and hence we can replace their 
elementary unit eigenvectors by arbitrary unitary linear combinations thereof. For example, 
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the four Lindblad operators which describe the T\ relaxation of the first spin may be replaced 
by the Hermitian operators: 

= v / OÏ532Í(|2)(0| + |0)(2| + |3)(l| + |l)(3|) = >/0l532Í^ 

Lg = V0.1532 | (|2)(0| - |0)(2| + |3)(1| - |1)(3|) = V0.1532|^ (25) 

L^? = v / 0l532|(|2)(0| + |0)(2|-|3)(l|-|l)(3|) = V0A532 ± 

L% = v/01532 1 (|2)(0|-|0)(2|-|3)(1| + |1)(3|) = vCÏ532Í^ z 2 

In a similar fashion, we may take those describing 7\ relaxation of the second spin to be 

L% = V0.1528±c^, Lj! = v/01528 i <j 2 , 

1 2 1 2_y (26) 

L£ = V0.1528i^, L£ = V0.1528|^ . 

In addition, the near-degeneracy of the (1,4) and (2,3) rates can be used to combine the 
associated Lindblad operators into four multiple-quantum T\ Lindblad operators based on 
the average rate: 

L£ = V0.0252|4^, = V0.0252±^ . 

By working through some examples, it may be seen that the sum of the Lindbladian 
superoperators for each of the three sets of four Lindblad operators above also causes all 
the off -diagonal entries of to decay with the same rate constant 1/Ti. This corresponds 
to nonadiabatic T 2 relaxation. Because does not act on its off-diagonal entries, it 
may also be seen that if one takes the Lindblads {(^R^^k)) 1 / 2 \k){k\ of the four diagonal 
entries of R^ and subtracts their superoperators from those for the off-diagonal, this must 
exactly cancel the nonadiabatic T 2 decay. Formally, however, the negative of a Lindbladian 
superoperator is not a Lindbladian superoperator, and in any case we do not really want 
to cancel the nonadiabatic T 2 , since it actually occurs. In order to avoid accounting for the 
nonadiabatic T 2 twice, it is nevertheless necessary to write down a set of Lindblad operators 
for it alone, without any T\ relaxation. Once again, on using the near equality of the 
diagonal entries of R^ to replace them by their average and taking suitable unitary linear 
combinations of the diagonal Lindblads \j){j\, we obtain 

L5f = V0.3312ÍJ, l^ 1 = VO.3312!^ 1 , 

2 ~ 1 2 (28) 

LT T f = v^33Ï2|a z 2 , L5f = v^33Ï2|a z V z 2 . 
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The Lindbladian superoperator of the first of these is obviously L^(g_ A ) = 0, and so need 
not be considered further. 

We now turn our attention to the diagonal entries of the 16 x 16 Zeeman relaxation 
superoperator 7?. zee , which we shall arrange in a 4 x 4 matrix of relaxation rate constants of 
the corresponding entries of the density matrix £^ e . It can be shown that this 4x4 matrix 



pzee 
Aldiag 



(00|^ zee |00) (01|ft zee |01) (02|7£ zee |02) (03|ft zee |03) 

(10|ft zee |10) (ll|^ zee |ll) (12|7£ zee |12) (13|7£ zee |13) 

(20|7?. zee |20) (21|7?. zee |21) (22|^. zee |22) (23|^. zee |23) 

_(30|ft zee |30) (31|ft zee |31) (32|ft zee |32) (33|^ zee |33)_ 



(29) 



is also a symmetric submatrix of the Choi matrix Ç associated with ]Z zee (up to sign), 
specifically 



ildiag 



(00|Ç|00) (00|Ç|11) (00|Ç|22) (00|Ç|33) 

(11|Ç|00) (11|Ç|11) (H|Ç|22) (11|Ç|33) 

(22|Ç|00) <22| C [11> (22|Ç|22) (22|Ç|33) 

(33|Ç|00) (33|Ç|11) (33|Ç|22> (33|Ç|33) 



(30) 



The diagonal of Rj^ g is thus the same as the diagonal of R^- Since we have already found 
a set of Lindblad operators which fully account for the effects of on £ z £ e , we will now 
focus our attention on its off-diagonal entries by defming a new matrix Rï^~ which is the 
same as R^ K save for its four diagonal entreis, which are set to zero as indicated below: 



pzee p zee 

£±T 2 — ildiag 



D\^((kk[R zee \kk) |fc = 3). 



(31) 



As implied by our notation, R^ contains all the information regarding T2 relaxation 
processes that is contained in our full, but diagonal, relaxation superoperator 7?. zee , and 
in a considerably more compact and easily understood form. Unlike R^, which acts on 
the column vector diag(/? z ^ e ) of diagonal entries by matrix multiplication, R^ acts on £_ z £ e 
simply by taking the products of all corresponding pairs of entries, one from each matrix, just 
as these entries are multiplied together in the full matrix- vector product 7?. zee |£a 6 ). This 
"entrywise" matrix multiplication, commonly known as the Hadamard product ||15|| , has 
already been shown to be a powerful means of describing "simple" T2 relaxation processes 
(that is, processes not involving cross-relaxation). The Hadamard product will be 
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denoted in the following as: 



àT = -m®aT = -[(j\pT\k)(j\R z %\k)] 3 . k=0 . (32) 



j,k=o ■ 



Another important property of the matrix Rj^ 2 is that, since the overall projected Choi 
matrix S_CS_ must be positive semidefinite, the same is true of the projection of —R^, anc ^ 
the 4x4 block of the 16 x 16 projection matrix £_ that acts on —BJr^ ^ s 



K = L- |II T = \ 



3 -1 -1 -1 
-1 3 -1 -1 
-1 -1 3 -1 
-1 -1 -1 3 



(33) 



where 1 denotes a column vector of four l's. Moreover, the Lindblad operators for T 2 
relaxation may be extracted directly from —E_R^ 2 E_ without direct reference to the full 
superoperator's projected Choi matrix. 

To see how this can be done, we first observe that given any two diagonal matrices A & 
C_ (assumed here to equal in dimension) with column vectors of diagonal entries a = diag (A) 
& ç = diag (C), respectively, we have 

ARC = (ill T Ç)0B = (aç T )0fi (34) 

for any other (not necessarily diagonal) matrix B_ of equal dimension. It follows that the 
action on a density operator p of any Lindblad operator L, with respect to a basis wherein 
its matrix L is real and diagonal, can be expressed in terms of Hadamard products as 

Up) = LqL — \ Uq_ — \ q_U = (L11 T L- ÍL 2 11 T - ill T L 2 )©£ = -R L &£, (35) 

where R L is called a Hadamard relaxation matrix for L. If múltiple Lindblad operators act 
simultaneously, the net Hadamard relaxation matrix is of course the sum of those associated 
with the individual Lindblads. 

Next, let us use the Lindblad operators for nonadiabatic T 2 relaxation given in Eq. (|28|) 
above to illustrate how we can go the other way, that is derive these Lindblads from the 
corresponding Hadamard relaxation matrix. If we let £™ = diag (L^) be the column 
vectors formed from the real diagonal entries of these Lindblad matrices and observe that 
their Hadamard squares £™ £™ = 0.3312 1 1 for i = 1, 2, 3, then this may be written as 

_pna _ «nal /«nal^T , p na2 / p nal \ T , p nal In nal\ T 

- \ {Lt? © It" 1 + írf © It? + ÍtT © ÍtT) I T 

- \ I {t™ 1 © It? + It? © It? + & & T 
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0.3312 
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(36) 



0.3312(E - §11 T ) = -0.3312 



111 
10 11 
110 1 
1110 



Noting that = E_ and that the row and column sums of E are zero, it may be seen 
that the projection —E_R£?E_ simply removes the last two terms involving the Hadamard 
squares from the above, which are proportional to 11 , leaving only —E_R^?E_ = 0.3312 i? 
behind. Because the vectors ££? % are mutually orthogonal and all their squared norms are 
Wi^ II 2 = 0.3312, upon normalization they actually become the eigenvectors of E_ that are 
associated with its one nonzero, but triply degenerate, eigenvalue of unity. 

From this we see that the nonzero entries of the diagonal Lindblad matrices are the 
entries of the eigenvectors of —E_R™ E_ times the square roots of their respective eigenvalues, 
much as general Lindblad matrices may be obtained from the eigenvalues and eigenvectors 
of the projected Choi matrix £Ç£_. The numerical vàlues of the entries oí R^, as extracted 
from the experimental superoperator in FIG. 0, are: 



75 zee 
£±T 2 



0.7890 0.7757 1.2872 
0.7890 0.5033 0.7426 
0.7757 0.5033 0.7283 
1.2872 0.7426 0.7283 



(37) 



It is easily seen that RJr^i n ^ e BJr* \ niust be centrosymmetric, and if we likewise symmetrize 
and subtract the above nonadiabatic T 2 Hadamard relaxation matrix, we get 



p ad 

£±T 2 



pzee p na 

ÜT2 £±T 2 



0.4274 0.4279 0.9560 
0.4274 0.1721 0.4279 
0.4279 0.1721 0.4274 
0.9560 0.4279 0.4274 



(38) 



The nonzero eigenvalues and associated eigenvectors of —E_R^ 2 E_ are (to within 1%) 





" i " 




- i - 




" 0" 


^0.9560, ^ 





. -i . 


V f 0.2913, i 


-i 
-i 
. i . 


\ (°- 1721 > 75 


i 

-i 

. 0. 



(39) 



which correspond to a system of three Lindblad operators for the adiabatic T 2 relaxation, 
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namely 

L·t 2 1 = V^9560 ^(2Í+d), kt = VoTmi^-ay, 

(40) 

These correspond to totally correlated, totally anticorrelated, and pure single-quantum, i.e. 
dipolar T 2 relaxation, respectively |ÜJ. 
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